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ABSTRACT 

Context. Spectral analyses of hot, compact stars with non-local thermodynamical equilibrium (NLTE) model-atmosphere techniques 
allow the precise determination of photospheric parameters such as the effective temperature (r e ff), the surface gravity (log g), and the 
chemical composition. The derived photospheric metal abundances are crucial constraints for stellar evolutionary theory. 

Aims. Previous spectral analyses of the exciting star of the nebula A 35, BD-22°3467, were based on He+C+N+O+Si+Fe models 
only. For our analysis, we use state-of-the-art fully metal-line blanketed NLTE model atmospheres that consider opacities of 23 el- 
ements from hydrogen to nickel. We aim to identify all observed lines in the ultraviolet (UV) spectrum of BD-22°3467 and to 
determine the abundances of the respective species precisely. 

Methods. For the analysis of high-resolution and high signal-to-noise ratio (S/N) far-ultraviolet (FUSE) and UV (FIST/STIS) ob- 
servations, we combined stellar-atmosphere models and interstellar line-absorption models to fully reproduce the entire observed 
UV spectrum. 

Results. The best agreement with the UV observation of BD-22°3467 is achieved at T tS = 80+10 kK and log g = 7.2 ± 0.3. While 
T eS of previous analyses is verified, log g is significantly lower. We re-analyzed lines of silicon and iron (1/100 and about solar abun- 
dances, respectively) and for the first time in this star identified argon, chromium, manganese, cobalt, and nickel and determined 
abundances of 12, 70, 35, 150, and 5 times solar, respectively. Our results partially agree with predictions of diffusion models for 
DA-type white dwarfs. A combination of photospheric and interstellar line-absorption models reproduces more than 90% of the ob- 
served absorption features. The stellar mass is M » 0.48 M 0 . 

Conclusions. BD-22°3467 may not have been massive enough to ascend the asymptotic giant branch and may have evolved directly 
from the extended horizontal branch to the white dwarf state. This would explain why it is not surrounded by a planetary nebula. 
Flowever, the star, ionizes the ambient interstellar matter, mimicking a planetary nebula. 

Key words, stars: abundances - stars: atmospheres - stars: evolution - stars: individual: BD-22 3467 - white dwarfs - 
planetary nebulae: individual: A66 35 


1. Introduction 

Abell 35 (or 2 ooo = 12 h 53 m 32 s .79, d 2 ooo = -22° 52' 22" 5 5) was dis- 
covered by Abell (1955, A55 24) and classified as a planetary 
nebula (PN). Abell (1966) characterized it (A6635, henceforth 
A 35, PN G303. 6+40.0) as a homogeneous disk PN with an an- 
gular size between 636" and 938". Figure 1 shows an image of 
A 35 with very long exposure time. Its size is about 17' in east- 
west and 14' in north-south direction. The nebula shape appears 
to be quite atypical for a PN. The bow-shock structure is sur- 
rounded by a symmetric emission area that does not agree with 
the usual bi-polar or ellipsoidal morphologies. 

Grewing & Bianchi (1988) discovered a hot companion to 
the visible nucleus in the optical (SAO 181201, Jacoby 1981). 
They suggested that it is an extremely hot DAO-type 1 white 

* Based on observations with the NASA/ESA Flubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26666. 

** Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. 

*** Tables 1 and 2, and Fig. 2 are available in electronic form at 
http : / /www . aanda . or g 

1 Spectral classification: characteristic H I and He II absorption lines. 


dwarf (WD). A35, BD-22°3467 (LWHya) is therefore a re- 
solved binary system (Jacoby 1981; De Marco 2009), composed 
of a WD star of spectral type DAO and a G8 III-IV-type compan- 
ion (Unruh et al. 2001; Strassmeier 2009). The latter dominates 
the spectrum at A <; 2800 A (Grewing & Bianchi 1988). 

Borkowski et al. (1990) explained the PN asymmetry (Fig. 1) 
as the result of an interaction between the PN and the surround- 
ing interstellar medium (ISM), the observed bow shock (the star 
is moving with v = 125 kms -1 ) being located at the equilibrium 
sphere of ram pressure of the accelerated PN gas and the stellar 
wind ram pressure. 

A mean (HlPPARCOS, Tycho-2, and UCAC-2 values) proper 
motion of ^ i a = -60.95 ± 1.5mas/yr and fis = -14.63 ± 
1.4mas/yr was measured by Kerber et al. (2008) in agree- 
ment with HlPPARCOS values (HIP 62905, Van Leeuwen 2007) 
of n a = -60.91 ± 1.62mas/yr^ = -13.41 ± 1.24mas/yr. With 
its high proper motion, BD-22°3467 passes through the com- 
plete visible nebula within about 16000 years (Fig. 1). Recently, 
Frew & Parker (2010) claimed A 35 to be a “bow shock nebula in 
a photoionized Stromgren sphere in the ambient ISM”. This sce- 
nario does not include a PN. Weidmann & Gamen (2011), how- 
ever, classified A 35 as a binary PN (? + G8 IV) with a bc-CSPN 
(corresponding to binarity for the cool CSPN). 


Article published by EDP Sciences 
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Fig. 1. A 35: narrow-band composite of Ha, TO in], and [Sll] (in to- 
tal 36h exposure time) by Dean Salman (Sh 2—3 13, http://www. 
sharplesscatalog.com/). N is up, E is left, the FOV is 20' x 20'. 
The (red) arrow indicates the projected motion of BD-22°3467 during 
the last 10 000 years. 

HlPPARCOS measured a parallax of 7.48 ± 1.55 mas 
(Perryman et al. 1997, D = 134^^ pc). This was corrected later 
by Gatti et al. (1998) to D = 163^ pc (this value was adopted 
by Herald & Bianchi 2002). Gatti et al. (1998) determined a 
separation of the hot and cool component of 13-28 AU at this 
distance. However, they adopted 160 pc as a minimum distance 
because previous estimates (see their Table 2) gave much higher 
values (e.g. 360 ± 80 pc from photometry measured by Jacoby 
1981). Van Leeuwen (2007) presented a validation of the new 
HlPPARCOS reduction. The improved parallax (8.38 ± 1.57 mas) 
is slightly larger with a similar error range. The corresponding 
distance of D = 1 19 42 * pc is now even smaller. 

Herald & Bianchi (2002) determined the atmospheric 
parameters of both components and found T e ([ = 80 kK, 
log g = 7.7 (cm/s 2 ) for the WD, and T e ff = 5 kK, log <7 = 3.5 
for its companion star. Their hot-component analysis of ul- 
traviolet (UV) spectra (FUSE 2 , HST/STIS 3 , and IUE 4 ) cover- 
ing 905 A $ 1 $ 3280 A was performed with TLUSTY and 
SYNSPEC (Hubeny 1988; Hubeny & Lanz 1992, 1995; Hubeny 
et al. 1994). They determined abundances for He, C, N, O, Si, 
and Fe (Table 4). For the analysis of the cool companion they 
used local thermodynamic equilibrium (LTE) models (Kurucz 
1991). 

State-of-the-art non-local thermodynamic equilibrium 
(NLTE) model atmospheres are fully metal-line blanketed and 
can consider all elements from hydrogen to nickel (Rauch 2003; 
Rauch et al. 2007). Therefore, we decided to re-analyze the 
available UV spectra of BD-22°3467 to determine abundances 
of hitherto neglected metals. The resulting abundance pattern 
can give clues to the evolutionary history of the DAO star. We 

2 Far Ultraviolet Spectroscopic Explorer. 

3 Hubble Space Telescope/Space Telescope Imaging Spectrograph. 

4 International Ultraviolet Explorer. 



Fig. 3. UV-optical-IR flux distribution of the BD-22°3467 binary sys- 
tem, fitted by atmosphere models for both components. 


begin with a brief description of the observed UV spectra and 
preparatory work (Sect. 2), followed by an introduction to our 
model atmospheres (Sect. 3). Then we describe our spectral 
analysis (Sect. 4) in detail and discuss the results in Sect. 5. 

2. Observations and interstellar absorption 

Hot post-AGB (asymptotic giant branch) stars have their flux 
maximum in the UV wavelength range. Most of the exhibited 
metal lines are located there. Metal lines of successive ioniza- 
tion stages allow one to evaluate the ionization equilibrium of 
the respective species. This is a very sensitive indicator for 5" e u. 
The error ranges of T e s and abundance determinations can be 
reduced by the analysis of many spectral lines. 

The dominant ionization stages of the iron-group elements 
in the relevant T e g and log g regime are V and VI (Fig. 2). Many 
strategic lines of these ions are located in the UV wavelength 
range. We therefore retrieved FUSE, HST/STIS, and IUE obser- 
vations from the MAST 3 archive (Table 1). 

Our spectral analysis is based on the high-resolution FUSE 
and STIS observations. The four STIS spectra were co-added to 
improve the S/N. Herald & Bianchi (2002) analyzed the same 
FUSE spectrum but only one (04GT02010) of the STIS spectra. 
The low-resolution IUE 6 observations were used in addition to 
determine the interstellar reddening (Eg-v, Fig. 3) as well as to 
verify the flux calibration. The STIS observation had to be scaled 
by a factor of 1 .995 to match the flux level of the IUE observa- 
tion. This offset is possibly due to bad centering of the exciting 
star during the exposure (Herald & Bianchi 2002). 

2. 1. Interstellar absorption and reddening 

The general procedure of determining the interstellar redden- 
ing is to normalize a theoretical spectrum as far as possible in 
the low-energy range because it is not significant there. In the 
case of BD-22°3467, the companion star dominates the ob- 
servation at wavelengths longer than ~2800A. Therefore, we 
constructed a theoretical, composite spectrum (with appropriate 

5 http://archive.stsci.edu/ 

6 International Ultraviolet Explorer. 
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Fig. 4. STIS observation around H I Ly a compared with our final model 
(Table 4) at different n H i. The dashed line is the WD’s atmospheric flux. 


stellar radii, R 2 d /R 2 = 7.43 x 10 5 ) using our final model spec- 
trum (Tf.fi = 80 kK. log g = 7.2) and a Kurucz model 7 8 * (T e g = 
5000 K, log g = 3.5, solar abundances logZ = 0.0) for the hot 
and cool component, respectively. This spectrum was then nor- 
malized to the 2MASS K flux (Fig. 3). The best fit (match- 
ing all GALEX, HlPPARCOS, and 2MASS magnitudes and the 
UV spectra) was achieved for E B _y = 0.02 ± 0.02, using the red- 
dening law of Fitzpatrick (1999) with the standardly = 3.1. Our 
value agrees within error limits with those of Herald & Bianchi 
(2002, Eb-v = 0.04 ± 0.01), who used the Rayleigh- Jeans tail of 
the UV observations (1000 A A T, 2800 A) only. 

To determine the column density of interstellar neutral hy- 
drogen, h H i, models with different « H i were compared with Ly a. 
We found a best fit for n H i = 5.0 ± 1.5 x 10 20 cm -2 (Fig. 4), 
which agrees with the value found by Herald & Bianchi (2002). 

2.2. Radial velocity of the hot component 

The STIS observation shows numerous sharp, isolated photo- 
spheric lines that are suitable to determine y ral j. On average 
(2 1 lines of N v, O iv, O v, Ar v, Cr v, Cr vi, Fe v, Fe vi, and Ni v), 
we measured u ra d = - 14.4 ± 0.7 km s _1 . 


3. Model atmospheres and atomic data 

In this section, we briefly describe the atomic data and the pro- 
grams used for our analysis. Details about the calculation of 
the stellar atmosphere models can be found in Sect. 3.1. The 
usage of OWENS for the calculation of the ISM absorption is 
described in Sect. 3.2. 

3.1. The photospheric model for the hot component 

Our model atmospheres (plane-parallel, chemically homoge- 
neous, in hydrostatic and radiative equilibrium) are calculated 
with TMAP X , the Tubingen NLTE model atmosphere package 
(Werner et al. 2003). We calculated a small grid of models that 
include opacities of 23 elements from H-Ni (Table 2). 


7 ftp : //ftp . stsci . edu/cdbs/grid/k93models/kp0®/ 

8 http : //astro . uni - tuebingen . de/~TMAP 


H-Ar are represented by “classical” model atoms (Rauch 
1997) taken from TMAD 7 , the Tubingen model atom database. 
For Ca-Ni, a statistical approach was applied to handle the 
large number of atomic levels and line transitions. We employed 
IrOnIC (Rauch & Deetjen 2003; Rauch 2003) to consider thou- 
sands of levels and millions of lines provided by Kurucz (2009, 
and priv. comm.) and the Opacity Project (Seaton et al. 1994). 

Kurucz‘s line lists are divided into so-called LIN (measured 
and theoretical lines) and POS (only measured, “good” wave- 
lengths) lists. For the model-atmosphere calculations, the LIN 
lists were used, to consider the total opacity properly. To cal- 
culate the emerging spectrum, POS lists were used to identify 
iron-group lines. Figure 5 demonstrates the difference between 
LIN and POS line lists. As an example for all species and their 
ions serves Crv (within the shown 1470- 1500 A interval). At 
present, Kurucz provides 73 222 LIN and 249 POS lines of Cr v, 
459 LIN and 6 POS lines in this wavelength interval. We ad- 
justed the Cr abundance to [Cr] = +3.1 to reproduce the two 
strongest Crv POS lines, Crv AA 1482.76, 1489.71 A. At this 
abundance, most of the prominent Crv LIN lines appear dis- 
placed or simply too strong. An exception is Crv A 1481.66 A, 
which matches the observation well. It can not be excluded that 
the Crv A 1490.24 A LIN lines in the model is the observed 
1490.66 A line, for instance. However, reliable abundance de- 
terminations are only possible using unambiguously identified 
POS lines (Sect. 4.3). 

In the framework of the Virtual Observatory (VO 10 ), all spec- 
tral energy distributions (SEDs, A - FT) calculated from our 
model-atmosphere grid are available in VO-compliant form from 
the registered VO service TheoSSA 11 provided by the German 
Astrophysical Virtual Observatory (GAVO 12 ). 

3.2. The ISM line-absorption model 

In the FUSE wavelength range (Fig. 6), numerous interstellar H 2 
and H 1 absorption lines hamper the analysis of the photospheric 
spectrum. We applied the program OWENS to calculate ISM 
line-absorption models. OWENS allows one to model different 
clouds with distinct radial and turbulent velocities, temperatures, 
column densities, and chemical compositions. It fits Voigt pro- 
files using ax 2 minimization. For more information on OWENS, 
see e.g. Hebrard et al. (2002) or Hebrard & Moos (2003). 

Figure 6 shows a small section of the FUSE spectrum of 
BD-22°3467. Only the combination of the photospheric and 
ISM-model spectra reproduces most of its spectral features. 
Our ISM model includes lines of Hi, Hi (J = 0-9), HD, 
Cii-iv, Ni-in, Nv, Oi, Ovi, Mgn, Ain, Sin-iv, P 11 , Sn-m, 
Svi, Art-n, Mnn, Fen, and Nin. The radial velocity of the 
interstellar low-ionization (I— II) atomic gas and H 2 is found to 
be D ia d = 6 ± 2 km s _1 . 


4. Analysis 

In this section, we begin with a verification of those atmospheric 
parameters that were already determined by Herald & Bianchi 
(2002) which may deviate because these authors used TLUSTY 
and considered He, C, N, O, Si, and Fe, while we employed 
TMAP with far more elements (Table 2). This might result in 

9 http : //astro . uni - tuebingen . de/~TMAD 

10 http : //www .ivoa.net 

11 http://dc.g-vo.org/theossa 

12 http://www.g-vo.org 


A109, page 3 of 12 


A&A 548, A109 (2012) 



Fig. 5. Section of the STIS spectrum (black line) compared with our TMAP model. The SED in the upper panel was calculated using Kurucz’s 
POS data (the strongest lines are identified at top), the SED in the lower panel using LIN data. This section contains mainly Fe V-VI and Mn V-VI 
lines. The strongest Crv LIN lines are identified at the bottom. 
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Fig. 6. Top : section of the FUSE observation showing the necessity of a detailed modeling of the ISM. The dashed line is the final, pure photospheric 
model. The full (red) line shows a combined stellar and interstellar model. Blue marks at the bottom denote interstellar H 2 absorption lines. Bottom: 
comparison of (pure photospheric) TMAP models with log g = 7.2 (thick, black) and log g = 7 .7 (thin, red), normalized to the flux of the log g = 7.2 
model at 999 A for comparison. The full lines are calculated with the determined reddening of Ejj^ = 0.04 (Sect. 2.1), the dashed lines (shifted in 
log/j by -0.4 for clarity) with E^ v = 0. 


a different atmospheric structure. We begin with models with 
parameters of Herald & Bianchi (2002) and fine-tune them if 
necessary. 


4.1. Effective temperature and surface gravity 

Herald & Bianchi (2002) determined T e g = 80 kK and 
log g = 1 .7 . TMAP test calculations showed that their T e g value 
matches all ionization equilibria well, e.g. N iv/N v, O iv/O v, 
and Mnv/Mnvi (Sect. 4.3). A lower logg i = 7.2 yields a 


significantly better fit agree of the theoretical Hen A 1640.42 A 
line profile with the observation in the line core, which is not 
matched at \ogg = 1 .1 (Fig. 7). For the line-profile calcula- 
tions, we used the Stark line-broadening tables by Schoning & 
Butler (1989a,b). The line wings are too broad at loggi = 7.7 
and somewhat too weak at loggr = 7.2. The same is true for 
He ii A 1084.94 A (Fig. 7). The decrements of spectral line se- 
ries such as the Hi Balmer (Rauch et al. 1998), Hi Lyman 
(transitions n—n' = 1-3 to 1-8), and Hen 2 -n' series (2-5 
to 2-11) are strongly dependent on loggr. The fit to the outer 
line wings of the latter two series in the FUSE spectrum is also 
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Fig. 7. Spectra of our final model compared with the FUSE and STIS 
observations around Hen A 1084.94 A (top, dashed: pure photospheric, 
full line: photospheric + ISM line absorption) and Hen A 1640.42 A 
(bottom, dashed: logg = 7.7, full line: logg = 7.2). The identified lines 
are marked. Most of the very strong absorptions in the top panel stem 
from H 2 . 

very well reproduced at log g = 7.2. Figure 6 shows Hen 2-5 
to 2-14 and the respective Hi 1-3 to 1-7 blends. A comparison 
of their theoretical profiles at logg = 7.2 and logg = 7.7 with 
the observation shows that they are significantly too broad at 
logg = 7.7. This effect cannot be compensated for by a lower 
Ep-v (Fig. 6). The interstellar Hi absorption is not dominat- 
ing the wings of these lines. The observed “shoulders” between 
the H 1 Lyman lines and, thus, the Lyman decrement is well re- 
produced at logg = 7.2 and is definitely missed at logg = 7.7. 
This is best visble between Ly £ and Ly 5 because the contamina- 
tion of the FUSE observation with ISM lines is obviously weak, 
logg = 7.2 agrees with the value of Herald & Bianchi (2002) 
within realistic error limits of 0.3 dex. Their assumed statisti- 
cal error (log cj = 7.7^'Jg) seems to be too optimistic. We finally 
adopted T t s = 80 ± 10 kK and log g = 7.2 ± 0.3 for our analysis. 

4.2. Hydrogen and helium 

We adopted [He] = -0.3 ([X] denotes log [mass fraction/solar 
mass fraction] for species X) from Herald & Bianchi (2002) and 
achieved a reasonable fit to the He 11 lines (Fig. 7). 

4.3. Metal abundances 

In our TMAP model-atmosphere calculations, we considered 
opacities of all elements from H to Ni except for Li, Be, B, 



Fig. 8 . Civ AA 1548.20, 1550.77 A lines in the STIS observation com- 
pared with our final model. Blue, thin: photospheric lines, red, thick: 
photospheric + interstellar lines, black, thin: interstellar lines. 


Cl, and K. We convolved our synthetic spectra with Gaussians 
to simulate the respective instrument resolution. By comparison 
with the observed spectrum of BD-22°3467, we determined for 
the first time photospheric abundances of Ar, Cr, Mn, Co, and 
Ni. All identified lines of these species in the available obser- 
vations (Sect. 2) were evaluated in our abundance analysis. No 
lines of F, Ne, Na, Mg, Al, Ca, Sc, Ti, and V were identified. 
For these, only upper limits could be derived by test models 
where the respective lines in the model emerge (at the abundance 
limit) from the noise in the observation. A detailed comparison 
of models with different abundances of individual metals showed 
that the typical abundance error is ±0.3 dex. In the following, 
we briefly describe some measurements (ordered by increasing 
atomic weight). 

Carbon : lines of stellar origin are not detectable in the spectrum 
of BD-22°3467. Narrow interstellar C iv AA 1548.20, 1550.77 A 
lines are present (Fig. 8) that virtually “bracket” the photo- 
spheric lines = -14.4 kms -1 , Sect. 2.2). They have their 
origin in two different interstellar “clouds” with different radial 
velocities and column densities (cloud 1: u ra d = -5. 4 kms -1 , 
log(nciv/(g/cm 2 )) = 16.8; cloud 2: v m d = -30.4kms -1 , 
log(nciv/(g/crn 2 )) = 14.6). The radial velocity of cloud 1 agrees 
with the velocity of the low-ionization ISM gas (Sect. 3.2). 
Interestingly, Schneider et al. (1983) determined t) ral j = -6.6 ± 
3.8 kms -1 for the nebula. This even suggests a nebular origin of 
cloud 1 while cloud 2 with its significantly higher velocity lies 
on the line of sight toward us. 

It is worthwhile to note that no other interstellar two- 
component absorption-line feature is identified in the spectrum 
of BD-22°3467. 

However, a photospheric component, even at low abun- 
dance, worsens the agreement between model and observa- 
tion in the shoulder between the interstellar lines (Fig. 8). 
Since the strongest Civ lines in the FUSE wavelength range 
(T/1948, 1107, 1168 A) become too strong for [C] > -2.9, we 
determine this as an upper limit. This value is about four times 
lower than that found by Herald & Bianchi (2002). 

Nitrogen: Niv A 1718.55 A and Nv AA 1238.82, 1242.80A are 
prominent in the STIS observation. Our final model reproduces 
the strengths of these lines and the ionization equilibrium at 
[N] = -1.7 (Fig. 9) well. 
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Fig.9. Niv A 1718.55 A and Nv AA 1238.82, 1242.80 A lines in the Fig. 10. Oiv AA 1338.61, 1343.51 A (left and middle) and 
STIS observation compared with our final model. Ovdl371.29A (right) lines in the STIS observation compared 

with our final model. 


Herald & Bianchi (2002) determined a much lower pho- 
tospheric N abundance (<10 -3 times the solar value) because 
of strong line wings forming at higher abundances. The reason 
may be a different approximation by the quadratic Stark effect in 
TMAP and TLUSTY/SYNSPEC. 

For the line-absorption cross-sections TMAP calculates 


O'iyM 


yfne 1 fij 
m t c Av d 


H(a, v). 


( 1 ) 


e and m e are electron charge and mass, c is the light velocity, 
Av d the Doppler width, fij the oscillator strength for the transi- 
tion i — » j, and H the Voigt function with 


T 

47tAvd ’ 


(2) 


and 



Fig. 11. Si IV AA 1393.76, 1402.77 A in the STIS observation compared 
with our final model. 


- Av 
Av d ’ 


(3) 


where Av = v - vo with the line frequency vo. 

According to Cowley (1970, 1971), the line-broadening due 
to the quadratic Stark effect is approximately given by 


r S tark = 5.5 X 10“ 5 ^= 

Vt 


(«! 


UP)2 


Z + 1 


(4) 


Oxygen: is present in the spectra with the ionization stages IV 
and V. We can reproduce the ionization equilibrium and the 
line strengths. Figure 10 shows Oiv AA 1338.61, 1343.51 A and 
O v A 1371.29 A at [O] = -2.6. This agrees well with the obser- 
vation and the value of Herald & Bianchi (2002). 

Fluorine: was not identified in the observation. We modeled 
Fv AA 1082.31, 1088.38 A and F vi A 1139.49 A. They would be 
visible in the observation at abundances higher than solar. 


where n e is the electron density, T is the temperature, is the 
effective principal quantum number of the upper level, and z is 
the effective charge seen by the active electron. Radiative (T ra d) 
and collisional damping by electrons is considered so that 

r = Trad + Tstark- (5) 


Magnesium: lines are weak and almost fade in the noise, 
even the strongest, Mg iv A 1683.003 A. Figure 5 shows 
Mg iv A 1490.433 A, which is not detectable in the FUSE obser- 
vation at solar abundance. Thus, we can only determine an upper 
abundance limit of [Mg] = 0.0. 


TMAP uses fj of 7. 8580 x 10 -2 and 1.5716X 10 -1 , and classical 
damping constants T ra d of 1.4401 x 10 +9 and 1.4494 x 10 +9 for 
N v AA 1242.80, 1238.82 A, respectively. 

A comparison of calculated Tst values (Sahal-Brechot, priv. 
comm.) has shown that Fj^ ey /r|^ _Brechot = 0.33 is a good 
approximation for the Nv resonance doublet within 10000 < 
T/K < 200 000. We introduced this factor and adopted Tst = 
3.0 • Tg“* ley for our analysis. 


Silicon: exhibits only the Si iv AA 1393.76, 1402.77 A reso- 
nance doublet (Fig. 11), other Si lines are too weak. A weak ISM 
contamination is present but we can determine [Si] = —2.1. This 
agrees well with the upper limit of Herald & Bianchi (2002). 

Argon: is strongly enhanced ([Ar] = +1.1). Similar to the cen- 
tral star of the PN Sh 2-216, LSV +46°21 (cf. Rauch et al. 
2007), Arvi AA 1283.96, 1303.89, 1307.38 A are identified in 
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Fig. 12. Section of the STIS spectrum around Ar VI AA 1283.96, 
1303.89, 1307.38 A compared with our final, composite photospheric 
and ISM spectrum. 


the STIS observation (Fig. 12) and reproduced by our final 
model. Our models show a strong Arvn A 1063.55 A line but 
this is blended by interstellar H 2 absorption and therefore not 
suitable for a T^g determination via the ionization equilibrium. 

A problem are lines like Arv A 1371.92 A (Fig. 10) which 
are much stronger in our models than observed. Since many 
ionization equilibria of different species such as Mn v— VI 
(see below) are well-matched, a higher T e g that would 
shift the ionization balance from Arv to Arvi can be ex- 
cluded. We checked the possibility that the oscillator strengths 
of Ar v lines that we use in TMAP (from the Opacity 
Project, Seaton et al. 1994) are too strong. The Arv triplet 
3p 3 3 S°-3p 43 P (Arv AA 1341.57, 1350.39, 1371.92 A) is found 
in Opacity Project data (far off) at 1455.77 A (without fine- 
structure splitting) with / = 4.2333 x 10 -2 . From this value, 
we calculated /-values of 4.7037 x 10 -2 , 1.4111 x 10 -2 , 
and 2.3518 x 10 -2 for the three fine-structure components. A 
comparison with recent work of Tayal et al. (2009, / 1371.92 = 
2.257 x 10 -2 and 2.834 x 10 -2 ) shows good agreement with 
the Opacity Project. Thus, the reason for the disagreement of the 
Ar v lines remains unknown. 

Chromium: is about 70 times enhanced ([Cr] = +1.8). In 
Fig. 14), Cr vi AA 1417.66, 1455.28 A are shown that are well- 
reproduced by our final model. Crv AA 1482.76, 1489.71 A 
(Fig. 5) would require a much higher [Cr] = +3.1 (Sect. 3.1). 
Since T e g and log g are well-determined (Sect. 4.1) and, thus, a 
strong change of the ionization balance toward Cr v is impossi- 
ble within the error limits, we have to consider these two lines 
as blends and not suited for a reliable abundance determination. 

Manganese: exhibits several Mnv and vi lines in the STIS 
wavelength range, e.g. (the strongest) Mnv AA 1359.24, 
1382.88, 1432.84, 1440.31, 1443.31, 1452.88, 1457.46, 
1480.75, 1486.50 A and Mnvi AA 1272.44, 1285.10, 1333.87, 
1345.49, 1356.85, 1391.17, 1391.22, 1548.43A (Fig. 13). All 
these lines are reproduced at [Mn] = +1.5. Since many lines 
of two successive ionization stages are suitable for spectral 
analysis, the Mn V-VI ionization equilibrium is one of the 
prime indicators for T e g for BD-22°3467. r e ff = 80kK is 
well-matched at log g = 7.2 (Sect. 4.1). 

Iron: has many lines of Fev and vi in the UV spectrum. 
In contrast to Herald & Bianchi (2002), who determined 
a 0. 1-0.5 times solar Fe abundance, we find it approximately 


Table 3. Mass, luminosity, radius, distance, and height above the 
Galactic plane of BD-22 D 3467. 


Post-AGB evolution Post-EHB evolution 


M/M 0 

0.52 ± 0.05 

0.48 ± 0.05 

log L/L q 

1 50 +035 

1 -0.45 

1 oo +0 - 50 
u -0.80 

R/Rq 

0 030 +ofl13 

-0.009 

0 029 +0 012 

D/ pc 

376 -M2 

361/;|f 

HI pc 

2331 33 

774+ 51 
z -62 


solar (e.g. Fig. 5). There are, however, some Fe lines that are 
stronger in our model than observed (see some Fevi lines in 
Fig. 14, bottom panel). Within the error range, we find the same 
Fe abundance as Herald & Bianchi (2002). 

Cobalt: is strongly enhanced. Co v AA 1331.50, 1349. 60A indi- 
cate [Co] = +2.2 (Fig. 14). 

Nickel: is also overabundant. Numerous Niv and vi lines are 
found in the UV spectrum. They are all reproduced at [Ni] = 
+0.7 (e.g. Figs. 5, 14). 

4.4. Mass, luminosity, radius, and distance 

We determined the stellar mass and luminosity from a com- 
parison with evolutionary tracks for hydrogen-rich post-AGB 
stars (Fig. 15, Miller Bertolami, priv. comm.) and for post-EHB 
(extended horizontal branch) stars from Dorman et al. (1993, 
Fig. 17). The results are summarized in Table 3. 

The stellar radius follows from R = \jG M/g, where G is the 
gravitational constant. We determined the spectroscopic distance 
D = 36U[3 3 pc (Table 3) following Wassermann et al. (2010). It 
is derived front 

D/pc = 6.6175 x 10“ 6 V^v/Fobs ■ M/M 0 x 10- lo «9 (6) 

with the astrophysical flux F v at the surface of the star as 
predicted by the model spectrum, and F 0 bs is the observed, 
extinction-corrected flux (we used F 1513 ^ to match the center of 
the GALEX far-ultraviolet band). D is about three times larger 
than the improved HlPPARCOS distance (D = I 1 9+ 2 * pc, Sect. 1) 
but it agrees with the photometric distance (360 ± 80 pc, lacoby 
1981) for the cool component. 

It is worthwhile to mention that there is no agreement with 
the HlPPARCOS distance at the maximum value of log <7 = 7.83 
given by Herald & Bianchi (2002) either. The calculated spec- 
troscopic distance (using M = 0.68 M e , R = 0.017 R 0 , and 
the respective model flux) is then D = 226 pc, about twice the 
HlPPARCOS distance. Although trigonometric parallaxes in gen- 
eral have to be regarded as hard constraints, orbital movement in 
binary systems such as BD-22°3467 may bear on the precision 
of HlPPARCOS measurements (cf. Rauch 2000, for LB 3459 = 
TYC9166-00716-1). Figures 6 and 7 clearly demonstrate that 
the spectra cannot be reproduced with a gravity of log 17 = 7.7 
or higher, while a good fit is achieved for log g = 7.2. The up- 
coming GAIA mission will provide a highly precise parallax for 
BD-22°3467 which will be able to rule out the uncertainties. 

5. Results and discussion 

We calculated a grid of TMAP NLTE models atmospheres 
for BD-22°3467, the exciting star of the ionized nebula A 35. 
We determined T e ff = 80±10kK and log <7 = 7.2 ± 0.3, and 
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Fig. 13. Strongest Mnv {top) and Mn VI lines ( bottom ) in the STIS wavelength range compared with our final model. 


thus confirm T e ff = 80 kK found by Herald & Bianchi (2002). 
Nevertheless, we found a somewhat lower log g than Herald 
& Bianchi (2002, log <7 = 7.7). For the first time in this object, 
abundances from detailed line-profile fits were determined for 
Ar, Cr, Mn, Co, and Ni (Table 4). For F, Ne, Na, Mg, Al, Ca, Sc, 
Ti, and V, upper limits could be derived. 

The element abundance pattern (Fig. 16) is probably the re- 
sult of the interplay of radiative levitation and gravitational set- 
tling. Compared to those objects in the DAO sample from Good 
et al. (2004, 2005) that have the same T e g and \ogg within er- 
ror limits, BD-22°3467 has a higher He and a lower Si abun- 
dance. Cr, Mn, Co, and Ni (not measured by Good et al. 2005) 
are strongly overabundant, likely due to the force of the radia- 
tion field. Radiative levitation may be a reason but the question 
remains why Fe is about solar. Detailed diffusion calculations 
that consider these elements are highly desirable investigating 
this phenomenon. A comparison with existing diffusion calcula- 
tions (Chayer et al. 1995a, b. Fig. 16) shows (within error limits) 
that the metal abundances partially agree with predictions for 
DA-type WDs. Only N, Mg, S, and Ca are more than one dex 
lower than predicted. 

Until recently, BD-22°3467 was believed to be the central 
star of a planetary nebula and, thus, a post-AGB evolution of an 
intermediate-mass star (0.8 M & Si Mj n i,j a i Si 8.0 M 0 ) appeared 
to be a matter of course. Frew & Parker (2010) showed that the 
visible, PN-like nebula may be a combination of shock excita- 
tion and photoionization of the ambient interstellar gas. 

BD-22°3467 provides enough high-energy pho- 
tons, (4.4 x 10 25 photons/s/cm 2 at energies >13.6eV 
and 6.4x 10 22 photons/s/cm 2 at energies >54eV calculated from 
our final model which give totals of A) 3 = 2.3 x 10 45 photons/s 
and A 54 = 3.3 x 10 42 photons/s, respectively) to ionize 

the surrounding interstellar gas (cf. Rauch et al. 2004). 
For A 54 /A 13 < 0.01, no Hen emission from the nebula 
can be expected. If we assume an electron temperature 
of T e = 10 000 K in the nebula, negligible extinction, d = 361 pc 
(Table 4), and an average angular diameter of 750", we derive 
that BD-22°3467 has a Stromgren sphere with a density 
of 16 cm -3 , log(F H / 5 /erg/s/cm 2 ) = -11.23, and an ionized 


Table 4. Results of this work compared with Herald & Bianchi (2002), 
whose He abundance we adopted. 



This work 

H&B 

NLTE code 

TMAP 

TLUSTY 

primary T e g- /kK 

80+ 10 

80 + 3 

log 9 

7.2 + 0.3 

-J 
-~o 
1 + 
op 

X H e(X G ) 

0.47 

0.3-0. 5 

Xc(X Q ) 

1.2 X 10 -3 

5.0 x 10 -3 -10 -2 

X N (X 0 ) 

1.8 X 10 -2 

<10 -3 

Xo(X Q ) 

2.5 X 10 -3 

5.0 x 10 -3 -10 -2 

X F (X Q ) 

SI 


X N e(X G ) 

<1 


X N a(X G ) 

<1 


X M g(X G ) 

<1 


Xai(X g ) 

<1 


X si (X G ) 

7.5 x 10 -3 

<10 -2 

Xp(Xo) 

<1.4 x 10 -3 


XsfXo) 

<1.2 x 10 -3 


XArtXo) 

12 


X C a(X G ) 

<1 


X S c(X G ) 

<1 


X Ti (X G ) 

<1 


Xv(X G ) 

<1 


Xcr(X G ) 

70 


XMn(Xo) 

35 


X F e(X G ) 

1 

0. 1-0.5 

X Co (X G ) 

150 


X Ni (X G ) 

5 


Eb-v 

0.02 + 0.02 

0.04 + 0.01 

log an i 

20.7 ±0.1 

20.9 + 0.1 

07 pc 

361 -137 

163 (Hipparcos) 

M a /M G 

0.48 + 0.05 

0.5 +0 ' 5 

RVRe 

0.029!°;^ 

0 0165 + ° 01 

u.wiuj_ 0006 

RVD*/R e / kpc 

0.080!»;“ 

0.1 + 0.01 

log L a /L 0 

1 oo +0 - 50 
1 ■ u -0.80 

1.0 +0.4 

secondary T t ff /kK 

5 

5 

log 9 

3.5 

3.5 


Notes. (a) Our values for post-EHB evolution. 


nebular mass of 0.5 M Q . The absolute H/l flux agrees very well 
with the value given by Acker et al. (1992, logF H ^ = -11.3). 
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1270 1271 1272 1273 1274 

x/A 


Fig. 14. Sections of the STIS observation containing (in addition to oth- 
ers) lines of Crvi (top), Co (middle), and Ni (bottom) compared with 
our final model. 

The nebular spectrum by Acker & Stenholm obtained for 
the catalog (Acker et al. 1992) is weak and, consequently, 
the Ffta/Fftp ratio is overestimated, and the derived extinction 
of c = 0.7 is very uncertain. F [OIII ] / 15007 ^/.Fh j s = 2 indicates 
low excitation. Both Nn AA 6548, 6584 A are visible and, with 
standard assumptions for the nebula but without any correc- 
tion for unseen higher ionization stages of nitrogen, its neb- 
ular abundance may be estimated to be 12 + logN/H = 7.7. 
The Si ii /! 67 1 6 A/Si ii T 673 1 A line intensity ratio is at the low- 
density limit and we can conclude that the nebular electron den- 
sity is n e < 100 cm 3 . In addition, we may estimate that the 
nebular oxygen and sulfur abundances are 7.9 and 6.4, respec- 
tively. It is worthwhile to note that the ionized mass of 0.5 M 0 
is of about the value usually assumed for the average PN mass 
of 0.2-0. 3 M 0 . The reason is simply the fairly low luminosity 
of BD- 22° 3467. Such a low-luminosity ionizing star makes a 
Stromgren sphere that may well mimic a genuine PN. It is thus 
possible that similar objects have been misclassified as PNe. 



Fig. 15. Location of BD-22°3467 in the log T e g -log g diagram com- 
pared with other DAO-type WDs taken from Gianninas et al. (2010, 
all symbols). Three of these (green squares) are used to compare their 
metal abundances with those of BD-22°3467 (Fig. 16). The post-AGB 
evolutionary tracks were calculated by Miller Bertolami (priv. comm.) 
labeled with the stellar mass (in M 0 ). 



5 10 15 20 25 

atomic number 


Fig. 16. Top', photospheric abundances of BD-22°3467 compared with 
solar values (Asplund et al. 2009), and mean values from three objects 
of Good et al. (2005, blue triangles) that are located within the error 
ellipse (cf. Fig. 15, green squares there). Arrows indicate upper lim- 
its. Bottom: comparison of number ratios compared with predictions 
of diffusion calculations for hydrogen-rich (DA-, blue squares) and 
helium-rich (DO-type, green triangles) WDs (Chayer et al. 1995a, b) 
with T eS = 80 kK and logg = 7.2. 

Because the time since its departure from the AGB is rel- 
atively long for BD-22°3467 due to its low mass (Table 3), 
all matter ejected previously on the AGB (as well as the PN) 
may have dissipated into the ambient ISM. Thus, the result of 
Frew & Parker (2010) does not stringently exclude a post-AGB 
evolution. However, an alternative stellar evolution may be pos- 
sible. Bergeron et al. (1994) discussed that the progenitors of 
most DAO-type WDs are post-EHB stars (M ~ 0.48 M a ). 
Recently, Gianninas et al. (2010) presented a spectral analy- 
sis of 29 DAO- and 1 8 hot DA-type WDs and found that their 
DAOs (Fig. 15) are hotter and more massive that previously 
determined. Accordingly, they concluded that most of these 
(“a rather mixed bag of objects”) are products of post-AGB 
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Fig. 17. Location of BD-22°3467 in the r e (f-logg diagram com- 
pared with sd(0)B-type stars close to the EHB taken from Edelmann 
(blue, open squares 2003) and DAO-type WDs taken from Gianninas 
et al. (red diamonds and green, filled squares 2010, cf. Fig. 15). The 
post-EHB evolutionary tracks are from Dorman et al. (1993, Y = 
0.288 « 7 g ) and labeled with the stellar mass (in M 0 ). 

evolution. BD-22°3467 is located amongst the low-mass stars 
of their sample (Fig. 17) and consequently, a post-EHB evo- 
lution cannot be ruled out. In this scenario, the stars evolve 
directly from the EHB to the WD state (AGB -manque stars, 
Greggio & Renzini 1990). Bergeron et al. (1994) proposed that 
a weak mass-loss may keep enough helium in the line-forming 
regions (Unglaub & Bues 1998, 2000) to exhibit strong Hen 
absorption lines. Figure 17 shows BD-22°3467 in compari- 
son to post-EHB evolutionary tracks (Dorman et al. 1993). One 
can estimate that an M ~ 0.48 M e track (with solar helium 
content on the horizontal branch, HB) will match the posi- 
tion of BD-22°3467. (Dorman's tracks for masses higher than 
0.475 M 0 and solar helium content are not fully calculated down 
to lower T e ff and higher log g.) On its way from the EHB to its 
present position, i.e., towards higher T e g and higher log < 7 , grav- 
itational settling, radiative levitation, and a weak wind are inter- 
acting and, thus, the atmospheric abundances may deviate from 
the HB abundances. 
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Table 1. Log of the UV observations of BD-22°3467. Table 2. Continued for Ca-Ni. 


Instrument 

Obsld 


Start time 

Exp. time 

Ion 

Super levels 

Super lines 

Sample lines 






(UT) 

(s) 

Caiv 

6 

16 

20291 









FUSE 

P1330101000 

2000-05-20 20:27:37 

4416 

Cav 

6 

21 

141956 

STIS 


04GT02010 

1999-04-17 21:14:49 

2050 

Ca VI 

6 

19 

114545 



04GT02020 

1999-04-17 22:37:03 

2800 

Cavil 

6 

21 

71608 



04GT02030 

1999-04-18 00:16:10 

2740 

Ca VIII 

6 

20 

9124 



04GT02040 

1999-04-18 01:52:54 

2740 

Calx 

1 

0 

0 

IUE 

SWP44121LL 

1992-03-06 04:50:23 

1440 

Sc IV 

6 

20 

15 024 


LWP22518LL 

1992-03-06 05:26:54 

1440 

Sc v 

6 

21 

261 235 








Sc vi 


19 

237 271 








6 

Notes. The FUSE observation was performed with the LWRS (low- 

Sc vn 

6 

20 

176 143 

resolution aperture, resolving power R = 

A/AA » 20000), STIS with 

Sc vm 

6 

21 

91935 

grating E140M ( R x 45 800), and IUE in low-resolution mode (R » 

Sc ix 

1 

0 

0 

300). 







Tiiv 

6 

19 

1000 








Tiv 

6 

20 

26 654 

Table 2. Statistics of H - Ar model atoms used in our calculations. 

Ti vi 

6 

19 

95 448 








Ti VII 

6 

20 

230618 








Ti VIII 

6 

21 

1 82 699 


Levels 



Levels 


TilX 

i 

0 

0 

ion 

NLTE 

LTE 

lines 

ion 

NLTE ] 

LTE lines 

Viv 

6 

19 

37130 

Hi 

12 

4 

66 

Navi 

43 

10 130 

Vv 

6 

20 

2123 

Hii 

1 

0 

- 

Navn 

1 

0 0 

V vi 

6 

19 

35 251 

He I 

5 

98 

3 

Mg hi 

1 

34 0 

Vvn 

6 

19 

112 883 

Hen 

16 

16 

120 

Mg iv 

31 

0 93 

Vvni 

6 

20 

345 089 

Hein 

1 

0 

- 

Mg v 

15 

37 18 

Vix 

1 

0 

0 

Cm 

6 

61 

12 

Mg vi 

1 

0 0 

Criv 

6 

20 

234 170 

Civ 

54 

4 

295 

Al in 

1 

6 0 

Crv 

6 

20 

43 860 

Cv 

1 

0 

0 

Aliv 

15 

2 0 

Cr vi 

6 

20 

4406 

N hi 

1 

65 

0 

Al v 

1 

16 0 

CrVH 

6 

19 

37 070 

Niv 

16 

78 

30 

Al vi 

14 

24 16 

Crvm 

6 

20 

132221 

Nv 

54 

8 

297 

Al vn 

1 

0 0 

Crix 

1 

0 

0 

Nvi 

1 

0 

0 

Si m 

3 

31 1 





Om 

3 

69 

0 

Si iv 

16 

7 44 

Mn IV 

6 

20 

719 387 

Oiv 

18 

67 

39 

Siv 

25 

0 59 

Mn v 

6 

20 

285 376 

Ov 

90 

36 

610 

Si vi 

1 

0 0 

Mn VI 

6 

20 

70116 

Ovi 

54 

8 

291 

Pm 

1 

9 0 

Mn vn 

6 

20 

8 277 

Ovu 

1 

0 

0 

Piv 

15 

36 9 

Mn VIII 

6 

20 

37 168 

Fm 

1 

6 

0 

Pv 

18 

7 12 

Mnix 

1 

0 

0 

Fiv 

1 

10 

0 

Pvi 

1 

0 0 

Feiv 

6 

20 

3102 371 

Fv 

15 

91 

31 

Siv 

6 

94 4 

Fev 

6 

20 

3 266 247 

Fvi 

12 

115 

16 

Sv 

21 

89 37 

Fevi 

6 

20 

991935 

Fvn 

1 

0 

0 

Svi 

18 

19 48 

FeVH 

6 

20 

200455 

Nen 

1 

33 

0 

S vn 

1 

0 0 

Fevm 

6 

18 

19 587 

Nem 

3 

43 

0 

Ariv 

1 

349 0 

Feix 

1 

0 

0 

Neiv 

3 

37 

0 

Arv 

32 

329 38 

Co IV 

6 

20 

552916 

Nev 

20 

74 

35 

Ar vi 

16 

168 21 

Co v 

6 

20 

1469 717 

Nevi 

1 

0 

0 

Arvil 

40 

1 12 130 

Co vi 

6 

18 

898484 

Nam 

1 

186 

0 

Ar VIII 

1 

0 0 

Co vn 

6 

19 

492913 

Naiv 

I 

237 

0 




Co vm 

6 

20 

88 548 

Nav 

8 

42 

9 

total 

741 2837 2514 

Co ix 

1 

0 

0 


Niiv 

6 

20 

2512561 

Niv 

6 

20 

2766 664 

Ni vi 

6 

18 

7408 657 

Ni vn 

6 

18 

4195 381 

Ni vm 

6 

20 

1 473 122 

Niix 

1 

0 

0 

total 

279 

884 

33 219 636 


Notes. The so-called “sample lines” are combined to “super lines” by 
IRONIC. 


A109, page 11 of 12 


log ionization fraction log n, / (1/cm 3 ) T / kK 


A&A 548, A109 (2012) 



-5 0 -5 0 -5 0 

I VI log m / (g/cm 2 ) 


-5 0 -5 0 -5 0 

I J MBHMa. VI log m / (g/cm 2 ) 

II vil 

hi ...... .r. yin 7“ e „ = 80 kK log g = 7.2 (cm/s ) 

IV » ........ ix 

— «- V 


Fig. 2. Temperature and electron density stratification along with ionization fractions of all elements in our final model. 


A109, page 12 of 12 


